US005908176A 


United States Patent [i9] 

Gilyard 


[ii] Patent Number: 5,908,176 

[ 45 ] Date of Patent: Jun. 1, 1999 


[54] IN-FLIGHT ADAPTIVE PERFORMANCE 
OPTIMIZATION (APO) CONTROL USING 
REDUNDANT CONTROL EFFECTORS OF 
AN AIRCRAFT 

[75] Inventor: Glenn B. Gilyard, Palmdale, Calif. 

[73] Assignee: The United States of America as 

represented by the Administrator of 
the National Aeronautics and Space 
Administration, Washington, D.C. 

[21] Appl. No.: 08/783,639 


[22] Filed: Jan. 14, 1997 

[51] Int. Cl. 6 B64C 13/18 

[52] U.S. Cl 244/203; 244/213; 244/219; 

244/195; 244/130; 701/3 


[58] Field of Search 244/203, 213, 

244/215, 219, 76 R, 194, 195, 76 A, 221, 
225, 90 R, 130; 701/3, 4, 8, 10, 11, 14 

[56] References Cited 

U.S. PATENT DOCUMENTS 


4,159,088 6/1979 Cosley . 

4,326,253 4/1982 Cooper et at. . 

4,485,446 11/1984 Sassi . 

4,723,214 2/1988 Frei . 

4,741,503 5/1988 Anderson et at 244/203 

4,797,829 1/1989 Martorella et at 701/3 

4,899,284 2/1990 Lewis et at 701/3 

4,932,611 6/1990 Horikawa 244/203 

5,457,630 10/1995 Palmer 701/3 


5,457,634 10/1995 Chakravarty 701/3 

5,478,031 12/1995 Piche 244/195 


Primary Examiner — William Grant 
Attorney , Agent , or Firm — John H. Kusmiss 

[57] ABSTRACT 

Practical application of real-time (or near real-time) Adap- 
tive Performance Optimization (APO) is provided for a 
transport aircraft in steady climb, cruise, turn descent or 
other flight conditions based on measurements and calcula- 
tions of incremental drag from a forced response maneuver 
of one or more redundant control effectors defined as those 
in excess of the minimum set of control effectors required to 
maintain the steady flight condition in progress. The method 
comprises the steps of applying excitation in a raised-cosine 
form over an interval of from 100 to 500 sec. at the rate of 
1 to 10 sets/sec of excitation, and data for analysis is 
gathered in sets of measurements made during the excitation 
to calculate lift and drag coefficients C L and C D from two 
equations, one for each coefficient. A third equation is an 
expansion of C D as a function of parasitic drag, induced 
drag, Mach and altitude drag effects, and control effector 
drag, and assumes a quadratic variation of drag with posi- 
tions 8; of redundant control effectors i=l to n. The third 
equation is then solved for h iopt the optimal position of 
redundant control effector i, which is then used to set the 
control effector i for optimum performance during the 
remainder of said steady flight or until monitored flight 
conditions change by some predetermined amount as deter- 
mined automatically or a predetermined minimum flight 
time has elapsed. 

6 Claims, 7 Drawing Sheets 
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IN-FLIGHT ADAPTIVE PERFORMANCE 
OPTIMIZATION (APO) CONTROL USING 
REDUNDANT CONTROL EFFECTORS OF 
AN AIRCRAFT 

ORIGIN OF INVENTION 

The invention disclosed herein was made by an employee 
of the United States Government and may be manufactured 
and used by or for the Government for governmental pur- 
poses without the payment of any royalties thereon or 
therefor. 

TECHNICAL FIELD 

The invention relates to using any or all redundant control 
effectors of an aircraft (subsonic or supersonic transports 
and fighters) to optimize actual in-flight performance of any 
specific aircraft configuration (particularly transport aircraft) 
during various flight conditions of different mission seg- 
ments and while experiencing changes in flight conditions. 
The term “redundant control effectors” is defined as any 
control surfaces, such as modulation and vectoring control 
systems and control systems that are in addition to the very 
minimum effectors (surfaces and systems) required to make 
any specific aircraft configuration capable of flight during 
any specific flight condition. 

BACKGROUND OF THE INVENTION 

Aircraft efficiency is a critical factor for airline profitabil- 
ity. A one percent performance improvement (1% fuel use 
reduction) for a fleet of transport aircraft can result in great 
savings, as much as approximately $100 million per year for 
the United States fleet of wide-body aircraft at current fuel 
costs (55 cents per gallon) and an additional $20 million per 
year for each ten cents per gallon increase in fuel price. 

Development of a significant amount of transport effi- 
ciency technology that started in the 1970’s and 1980’s has 
continued into the 1990’s. An aircraft energy efficiency 
(ACEE) program explored maneuver load control, elastic 
mode suppression, gust load alleviations (J. F. Johnson, 
“ Accelerated Development and Flight Evaluation of Active 
Controls Concepts for Subsonic Transport Aircraft — Volume 
I—LoadAlleviationlExtended Span Development and Flight 
Tests ,” NASA CR- 172277, 1984 ), relaxed static stability 
(W. A. Guinn, J. J. Rising and W. J. Davis, “ Development of 
an Advanced Pitch Active Control System for a Wide Body 
Jet Aircraft ,” NASA CR-172277, 1984) and reduced- area 
horizontal tail design (J. J. Rising, “ Development of a 
Reduced Area Horizontal Tail for a Wide Body Jet Aircraft," 
NASA CR- 172278, 1984). An advanced F-lll fighter tech- 
nology integration Mission Adaptive Wing (MAW) program 
developed and demonstrated variable -camber control for 
optimization of cruise and maneuver flight conditions 
(‘ Advanced Fighter Technology Integration F - 111 Mission 
Adaptive Wing (ITAR), NASA CP-3055, 1990; E. L. Friend, 
“Fligh t Buffet Characteristics of a Smooth Variable-Camber 
Mission Adaptive Wing for Selected Wing Flap Deflections 
(ITAR), NASA TM-4455, 1993). 

A member of the Airbus Industries team performed pre- 
liminary design work and wind-tunnel testing for imple- 
menting variable camber as a redundant effector in the 
A330/A340 aircraft. Benefits of variable camber include 
improved aerodynamic efficiency through improved lift-to- 
drag ratio (L/D), increased Mach number (M) capability, 
improved buffet boundary, increased operational flexibility, 
reduced structural weight and reduced fuel burn. Variable 
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camber has also increased aircraft development potential. 
However, the proposed use of a variable camber design as a 
redundant effector did not include development of a real- 
time adaptive performance optimization methodology. 

5 J. J. Spillman has provided an excellent dissertation 
relative to the fundamentals of variable camber as applied to 
transports (“The Use of Variable Camber to Reduce Drag, 
Weight and Costs of Transport Aircraft,” Aeronautical 
Journal, Vol. 96, January 1992, pp. 1-9). 

10 American manufacturers have also been actively involved 
in efficiency enhancement and have explored and imple- 
mented fixed-point rerigging of redundant control effectors 
to minimize airframe drag (“Long Live the Leviathan,” 
Flight International, Sep. 15-21, 1993, pp. 30-31 and Guy 
Norris, “New MD-11 Update Revealed,” Flight 
15 International, Dec. 21, 1994-Jan. 3, 1995, page 9). 

The literature is replete with reports documenting trajec- 
tory optimization algorithms and their benefits relative to the 
economics of commercial transports. In fact, all large trans- 
ports currently being produced have computer-based Flight 
20 Management Systems (FMS) that optimize the aircraft tra- 
jectory to minimize cost as a function of flight time and fuel 
price. However, the common basis for these algorithms are 
models of performance -related aspects of the particular 
aircraft model under specific flight conditions. As a result, 
25 the optimal trajectory is only as good as the onboard FMS 
models. In addition to the baseline onboard model having 
less than perfect accuracy, airframe and propulsion system 
degradation are factors which affect model accuracy. 

NASA’s Dryden Flight Research Center (DFRC), at 
30 Edwards, Air Force Base, California is active in transition- 
ing performance improvement technology to transport air- 
craft (Glenn B. Gilyard and Martin D. Espana, “On the Use 
of Controls for Subsonic Transport Performance Improve- 
ment: Overview and Future Directions ,” NASA TM-4605, 
35 1 994). Realizable performance benefits are smaller for trans- 
port aircraft than for high performance fighter aircraft. The 
designs of most transports have already been highly refined 
for good performance under a steady-state cruise flight 
condition. An algorithm developed on a Performance- 
40 Seeking Control (PSC) program was useful as an early 
demonstration of the benefits to be accrued on performance 
aircraft with detailed models available, though not suitable 
for implementing further performance optimization on 
transports primarily because of the fact that the algorithm 
45 was heavily based on a priori model data and absolute 
measurement accuracy (Glenn B. Gilyard and John S. Orme, 
“Performance Seeking Control: Program Overview and 
Future Directions," NASA TM -4531, 1993). Consequently, 
DFRC is exploring the application of measurement-based 
50 APO for performance improvement on transports using 
redundant control effectors. 

Adaptive Control Background 

Application of adaptive control to aircraft has been ongo- 
ing for more than 30 years with varying degrees of success. 
55 These applications have often centered on handling quality- 
related control system improvements, which often involve 
optimizing a very subjective, often ill-defined criteria typi- 
cally involving handling qualities, for example, pilot ratings. 
Because of the subjective nature of handling qualities, 
60 adaptive control techniques are not necessarily well-suited 
approaches to the problem. Also note that in many flight 
control applications, use of adaptive techniques has led to 
safety concerns about gain and phase margin reductions. 
Such reductions have contributed to stability and control 
65 problems. 

As such, adaptive control, as applied to flight control, has 
not found wide acceptance within the aerospace community. 
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Lack of interest in adaptive control is partially caused by the 
satisfactory results that have been obtained using conven- 
tional design techniques and by lack of an overriding reason 
to obtain similar results by using more complex techniques. 

Application of adaptive control is particularly advanta- 
geous when the optimization objective is well defined and 
there are significant unknowns about the aircraft and its 
operation. Application of adaptive optimal control to quasi- 
steady performance optimization has clear benefits that are 
not achievable in control design processes that are tailored 
to handling qualities issues. Quasi-steady performance opti- 
mization has well-defined objectives (i.e., minimize drag). 
For this reason, adaptive optimal control is well-suited to 
performance optimization. In addition, application of adap- 
tive optimal control, using a measured performance metric, 
is insensitive to modeling inaccuracies and measurement 
biases. Because low frequency constrained maneuvers are 
proposed, stability- and control-related safety issues and 
affects on ride qualities are greatly minimized. 

Regarding work of the Airbus Industrie team and Ameri- 
cans involved in exploring variable-camber performance 
optimization, neither side has devoted serious attention to a 
real-time Adaptive Performance Optimization (APO) algo- 
rithm. The Americans used either predetermined deflection 
schedules or a real-time, trial- and-error approach for camber 
control {Advanced Fighter Technology Integration F- 111 
Mission Adaptive Wing (ITAR), NASA CP-3055, 1990, and 
P. W. Phillips, and S. b. Smith, AFTI/F-111 Mission Adaptive 
Wing (MAW) Automatic Flight Control System Modes Lift 
and Drag Characteristics, AFFTC-TR-89-03, May 1989). 
In the case of the Airbus Industrie team, only model-based 
or experimentally determined scheduling was briefly men- 
tioned as a means of camber control (J. Szodruch and R. 
Hilbig, “Variable Wing Camber for Transport Aircraft,” 
Progress in Aerospace Sciences, Vol. 25. No. 3, 1988, pp. 
297-328). 

Application of adaptive optimal techniques to perfor- 
mance optimization does not require accurate models or 
absolute measurements. The adaptive optimal approach is 
based on real-time estimation of gradients of performance 
measures to control variables. These gradients are based on 
flight measurements and not based on predictions, except as 
noted herein below. In addition, because gradients are used, 
the approach is insensitive to measurement biases. 

An adaptive optimal approach is ideally suited for use on 
operational “fleet” aircraft where there is uncertainty in the 
aircraft model and absolute measurement accuracy. 
Likewise, adaptive performance techniques have a valuable 
role for commercial aircraft where small benefits over a 
20-30 yr service life can produce significant cost savings. 

Many issues enter into the performance optimization 
problem for subsonic transport aircraft. Foremost, there 
must be the potential for optimization, which implies redun- 
dant control effector capability (i.e., more than one means of 
trimming out the forces and moments to maintain a steady- 
state flight condition). Most aircraft have significant capa- 
bility in this area although taking advantage of this capa- 
bility is a different issue. Performing optimization from a 
condition that is already fine-tuned (based on wind tunnel 
and flight testing) places increased demands on high-quality 
instrumentation to sense small differences in an unsteady 
environment. 

APO compensates for all unique characteristics of the 
aircraft flight control system by continuously feeding back 
measurements of parameters that reflect the optimization 
objective, such as fuel flow (minimize) or velocity 
(maximize). A specific example is the use of symmetric 
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aileron deflection in an algorithm applied to optimally 
recamber the wing for all aircraft configurations and flight 
conditions to optimize (maximize) lift-to-(drag ratio which 
can in turn be utilized to produce: minimum fuel flow, 
5 maximum Mach number, maximum altitude or maximum 
loiter time. 

A feasibility study by Glenn Gilyard explored a prototype 
adaptive control law on a high-fidelity, nonlinear simulation 
of a first-generation wide-body jet aircraft that optimized 
10 wing-aileron camber for minimum aircraft drag at a given 
flight condition (‘ Development of a Real-Time Transport 
Performance Optimization Methodology,” NASATM-4730, 
1996). This technology applies tc selected current genera- 
tion aircraft and could be a requirement for future designs, 
15 such as proposed new large aircraft. 

Challenge for Present and Future Aircraft 

The APO of the present invention could play an important 
role in improving economic factors for the operation of 
aircraft. The challenge to in-flight performance optimization 
20 for subsonic and supersonic transport aircraft and fighters is 
the identification and minimization of very low levels of 
incremental drag. This is the key technological challenge 
because in order to provide an effective performance opti- 
mization algorithm, identification of incremental drag levels 
25 of one percent or less are required. 

STATEMENT OF THE INVENTION 

The present invention addresses the practical application 
of real time (or near real time) in-flight Adaptive Perfor- 
30 mance Optimization (APO) using a multiengine transport as 
an example and concentrates on onboard measurement and 
calculation of incremental drag from forced response 
maneuvers of redundant control effectors. 

The method is based on using an Inertial Navigation 
System (INS) which may include blending with Global 
Positioning System (GPS) information that can produce very 
accurate linear and angular displacement, velocity and 
acceleration measurements. Along with other more conven- 
tional measurements, such as true airspeed available from an 
40 Air Data Computer (ADC), some of this INS and GPS 
information is used to calculate winds and angle of attack, 
a. Thrust, T,-, is estimated from a prepared map in the form 
of a table made from a representative engine as a function of 

M, P c and EPR- where, M is Mach number obtained from the 

45 

ADC, P s is static pressure obtained from the ADC, and EPR 
is engine pressure ratio obtained from engine measurement. 
In other words, thrust is estimated from a representative 
engine model as a function of calculated variables (M,P) and 
engine measurement values of EPR. While a thrust map may 
not be particularly accurate in absolute terms, it is typically 
very accurate in relative terms, i.e., it accurately reflects 
small changes about the trim (nominal) throttle position. 
The following lift and drag equations are then used to 
55 calculate lift and drag coefficients, C L and C^,: 

Q=Lift/(^)={WA z//7 -2J, sin (a-ViMqS) (1) 

C D =Drag/(qS)={I.T i cos (a-r i)-WA xfp }/(qS) (2) 

where: q=M 2 P s dynamic pressure, lb/ft 2 from ADC, 
S=wing aircraft reference area, ft 2 , known constant, 
a= angle of attack, radians, from aircraft measurement or 
calculation, 

W=aircraft gross weight, is calculated based on takeoff 
65 weight and fuel flow, 

A zf = acceleration normal to flight path, g, (positive 
upward g’s) from INS, 
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A^f =acceleration along the flight path (positive forward 
g’s) from INS, 

T|= inclination ( known constant) of engine thrust relative 
to aircraft fuselage, radians and 
Z=summation of the bracketed quantity over the number 
of engines of the aircraft from i=l to N. Thrust related 
ram effects are assumed aligned with the gross thrust 
axis and are included in the other terms. 

Next calculate the following unknown variables: 

Cdo, Koi CijninCj^Cj)^ + , K t s and S l0 pt s 

in the following equation 

C D = (3) 

C Do + K 0 \C l - CfjninCj) } 2 + C Dm AM + C[) Hp A//p + ZKi(di - 6 iopt f 

where CD z?0 =minimum drag coefficient, 

C Lminc= C L at minimum CD, 

C z?w =coefificient of drag caused by Mach number, 

Co Hp = coefficient of drag caused by altitude, 

AM=change in Mach number from the time of turning on 
the APO, 

AHp=Hp change in altitude from the time of turning the 
APO on, 

K 0 and K, are drag equation coefficients, 
deposition of a redundant control effector i, where i is an 
integer from 1,2 ... n, 

8; o/ , r = 0 ptimal position of redundant control effector i, 
Z=summation over the number n of redundant control 
effectors in use from i=l to n. 

Equation (3) is a mathematical expression for C D that is a 
function of aircraft parasitic drag, induced drag, Mach and 
altitude drag effects, and aileron drag. That equation 
assumes a quadratic variation of C D with 8,-, and 8 ( - can 
include one or more redundant control effectors of the 
aircraft, such as symmetric ailerons, flaps, slats, spoilers, 
stabilizer, elevator, thrust, center-of-gravity and thrust 
vectoring, and propulsion system geometry and controls for 
longitudinal effectors and rudder, ailerons (normal or 
differential) and differential thrust lateral or directional 
effectors, so long as the designated effectors satisfy the 
definition of a redundant control effector for the condition of 
the specific flight segment to be optimized, such as but not 
limited to ascent, cruise, and descent. 

If tables for estimating thrust as a function of EPR are not 
available, equations (1) and (2) can be expanded as follows: 

T=T 0 +K t (EPR) @Mach & altitude«constant 

T a = to be identified 

Kj=d(T)/d(EPR)=constant to be identified 

C L =Lift/(^)={(U0A z/p -2(r o+ ^Pi?)) Sin (a-vi)}/(qS) 

C D =Dmgl (<yS)={2 (T 0 +K T (EPR)) cos (a -ri)-(W)A xfp }/(qS) 

C L and C D are now algebraic expressions which are 
included in the expansion of the expression for in 
equation (3) as follows: 


{E(7-„ + K T (EPR))cos(a - r,) - ( W)Aj p }/(qS) = 

Cdo + - Wo + K t (EPR)) sin(a - q)}/(qS) - C tai „c D } 2 + 

5 C Dm AM + Co H ^ AH P + T.Ki(6i - S^,) 2 

If the engines behave similarly, only two additional 
unknowns (T 0 and K r ) are added to the solution. To the 
io degree engines behave differently, additional pairs of T 0 and 
K t are added as unknowns to the solution. 

If a priori values of K 0 and C LminC are available, they can 
be utilized to reduce the number of°unknowns to be deter- 
mined. The algorithm may run either continuously in search 
15 of any change of optimal effector position h iopt or until a 
criteria has been met to assure optimal flight conditions have 
been reached (e.g., the optimal control position remains 
constant) and a new sequence of iterations may be called 
only at times when a criteria of change in any one or more 
20 flight variables has been met selected from a group com- 
prising incremental change in speed, altitude, or elapsed 
time. The flight crew may also call for a new sequence of 
iterations at any time, but should not do so unless there is 
some clear indication, such as from the instrument panel, 
25 that an optimal flight condition has changed perceptibly. 
Otherwise, since the forced response of the redundant effec- 
tors may themselves induce drag while the optimum h iopt is 
being sought, it is best for the crew not to make a judgement 
except where it is clear a deviation from the optimum has 
30 occurred. In any case, the forced-response maneuvers should 
be small with a raised-cosine excitation form of 100 to 500 
seconds to collect 100 to 500 or so sets of data at one set or 
so per second over the total excitation time at a rate of one 
to ten sets per second or so. (The rate can range from 1 to 
35 10 sets per second.) This form for the excitation so satisfies 
the smoothness criteria as to be imperceptible to flight crew 
and passengers and allows for collecting a sufficiently high 
number of data sets required for the accuracy of optimiza- 
tion estimates desired. 

40 The novel features that are considered characteristic of 
this invention are set forth with particularity in the appended 
claims. The invention will best be understood from the 
following description when read in connection with the 
accompanying drawings. 

45 BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. la illustrates a multiengine wide-body aircraft hav- 
ing many control effectors, a number n of which may be 
employed as redundant control effectors as the term is 
5Q defined, depending upon what flight condition in steady state 
has been established. 

FIG. lb illustrates how the Adaptive Performance Opti- 
mization (APO) feature of this invention is functionally 
applied to one redundant effector in one example. 

55 FIG. 2 illustrates a raised-cosine excitation form preferred 
to be used for forced-response maneuvers of selected redun- 
dant control effectors in the APO methodology of the present 
invention. 

FIG. 3 is a graph of variation of thrust required for 
60 trimmed flight as a function of symmetric aileron displace- 
ment and shows a thrust of 24.03x10 3 lb, an improvement 
of 364 lb (1.5%) at approximately 4.5° symmetric aileron 
deflections. 

FIG. 4 is a flow chart on an operational implementation of 
65 the APO algorithm. 

FIG. 5 presents a time history of forced excitation 
response cycle with both altitude- and velocity -hold turns on 
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during drag minimization maneuver using symmetric aile- 
rons as the sole redundant effector. 

FIG. 6 presents variation of calculated drag coefficient 
(with no corrections) as a function of symmetric outboard 
aileron which does not present a clear picture of the mini- 
mum drag point. 

FIG. 7 presents change in drag coefficient (corrected for 
C L changes) as a function of symmetric outboard aileron 
command, deg. 

FIG. 8 presents change in drag coefficient (corrected for 
C L and Mach number changes) as a function of symmetric 
outboard aileron command, deg. 

FIG. 9 presents change in drag coefficient (corrected for 
C L and Mach number changes) as a function of symmetric 
outboard aileron command, deg (with resolution on param- 
eters used in the analysis). 

FIG. 10 presents variation of total thrust with throttle 
lever for the data of the time history of FIG. 5. 

DETAILED DESCRIPTION OF THE 
INVENTION 

An optimal control effector combination can only be 
determined when there is an excess of control effectors 
defined hereinbefore as redundant control effectors b i where 
i is an integer from 1 to n corresponding to how many of the 
excess effectors are to be used in the APO methodology for 
a particular flight condition. For most aircraft at a level 
cruise condition, there are usually only two control effectors 
with no excess. For example, in a multiengine aircraft shown 
in FIG la, level flight at constant Mach and altitude 
(assuming the aircraft is trimmed up in the lateral and 
directional axes) requires only two control effectors, typi- 
cally thrust, T, (the throttle, not shown) and horizontal 
stabilizer 10. The throttle is used to counter balance the drag 
so that a given level of energy is maintained and the 
stabilizer 10 controls the balance between kinetic and poten- 
tial energy. 

The unique combination of throttle and horizontal stabi- 
lizer required to produce straight and level flight at constant 
Mach and altitude also results in a unique angle of attack, a,-. 
Many aircraft have elevators 11 but they are a minor force 
and moment generator relative to the horizontal stabilizer 
10. In many instances the elevators are coupled to the 
horizontal stabilizer motion. As such, the elevators are not 
an independent control effector. In such cases, the term 
stabilator refers to one control effector, namely the combi- 
nation of effectors 10 and 11. Other effectors shown are slats 
13, flaps 14, and rudder 15. 

Assuming that there is an additional control effector 
(more than the requred two discussed above) that can 
provide forces and moments in the aircraft’s longitudinal 
and pitch axes, the additional control effector is referred to 
as a redundant control effector since it can provide forces 
and moments similar to those of the throttle and horizontal 
stabilizer. For the sake of discussion, assume this redundant 
control effector comprises symmetric outboard ailerons 12 
on each wing which can be used to produce a change in wing 
camber and thus a wide range of wing aerodynamics. 

With this one redundant control effector, there are now a 
virtually infinite number of combinations (at least math- 
ematically speaking) of settings (positions) of the throttle, 
horizontal stabilizer, and symmetric ailerons which can 
produce a trimmed flight condition (level flight at constant 
Mach and altitude), each combination having a unique angle 
of attack. However, each of these combinations is not 
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equally good in the sense of providing minimum aircraft 
drag. Even though there can theoretically be multiple drag 
minima, it is more likely that there is only one pronounced 
minima over a range of small changes in the combinations 
5 which can produce a minimum drag configuration. A method 
for finding this one unique combination of redundant control 
effectors and resulting angle of attack a that the present 
invention addresses, namely an Adaptive Performance Opti- 
mization (APO) algorithm using at least one redundant 
10 control effector while in a steady flight condition, such as a 
climb, cruise or descent condition. 

Of course instead of there being only one additional 
redundant control effector there can be many, in which case 
a minimum drag for each may be sought by separate 
15 sequential excitations, each in accordance with the method- 
ology of this invention described herein with reference to 
equations (1), (2) and (3), but the set of minima, although 
good in the sense of providing a better “minimum” drag than 
without the methodology, would not guarantee a true opti- 
20 mal configuration. It would be preferable to integrate the 
excitations (forced response of redundant effectors) in the 
methodology using combinations of small perturbations 
affected simultaneously as stated in equation (3). 

The following describes the general sequence of events to 
25 obtain the optimal control effector combination, but first it 
should be noted that since the performance optimization 
method described is best suited to steady flight, the follow- 
ing autopilot and navigation modes should be engaged: 
altitude-hold, Mach- (or velocity-) hold, and heading-hold. 
30 In the event that the aircraft is not equipped for any of these 
modes, the pilot must perform the functions necessary for 
steady flight as best he can, although the latter case could 
result in less than true optimal b iopt being determined due to 
degraded quality in the forced response maneuver and 
35 analysis of the redundant control effector(s) resulting from 
less than perfect steady flight. The methodology will now be 
described. 

FIRST 

Forced response excitation must be provided for each of 
40 the redundant control effectors, i=l,2 . . . n. In the case of 
only one redundant control effector i=l, as in the example of 
the symmetric outboard ailerons 12 used as one redundant 
effector discussed above with reference to FIG. la, that 
effector must be excited (moved to a new position) using the 
45 smooth raised-cosine excitation form shown in FIG. 2 and in 
FIG. lb sufficiently for the resulting aircraft and control 
surface motion to provide sets of information (typically 100 
to 500 sets at the rate of one to ten sets per second for 
raised-cosine excitation of a time interval (period) equal to 
50 100 to 500 sec. from which the optimal position h iopt is to be 
determined. 

A passenger ride quality constraint is that the motions 
resulting from the forced response excitation be impercep- 
tible (less than~0.02 g) to the average passenger. This 
55 requires both a smooth and slow excitation such that the 
autopilot (or pilot) can respond to the excitation in an effort 
to control a steady altitude, Mach (velocity) and heading. 
The raised-cosine excitation form of FIG. 2 satisfies the 
smoothness requirements for excitation periods (intervals) 
60 of 100 to 500 seconds. An interval in that range should allow 
for ride quality constraints to be satisfied while collecting 
data sets required by the methodology of equations (1), (2) 
and (3) at rates of from one to ten sets per second. The 
raised-cosine excitation function of FIG. 2 is defined by the 
following equation: 

6~(A){1.0-cos ({2ot/P}TIME)}/2 (4) 
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where: A=maximum amplitude of the forced response 
(+or-) maneuver, 

P=period (duration) of the maneuver, 

TIME =starts at 0 and runs until the end of the period is 
reached, 
jt-3. 141592654 

Other excitation forms may, of course, be used to effectively 
obtain the optimal solution through analysis of all of the sets 
of data. However, the smooth raised-cosine form insures the 
smoothness that is required for passenger aircraft applica- 
tions and for all other aircraft as well. Thus, depending on 
the application, other excitation forms may be used. 

In any case, the excitation form for the maximum 
amplitude, A, of the forced response must be selected to 
cover the point of optimality, i.e., minimum drag which 
translates into minimum thrust required as a function of 
symmetric outboard ailerons at h iopt shown in FIG. 3 as 
equal to 4 to 5° for the example under discussion of 
optimizing cruise at 37,000 ft. That maximum excitation 
amplitude needed may be determined empirically. 

If there are multiple redundant control effectors, each 
must cycle through an excitation forced response maneuver. 
Analytically, the best means of determining the optimal for 
more than one redundant control effector (n>l) is to provide 
simultaneous excitation of the n redundant control effectors, 
as noted hereinbefore. However, the excitation of each of the 
redundant control effectors must be independent (i.e. 
different) in some way to insure identifiability of the separate 
control effectors. 

To achieve separability, the excitation signals used for 
each effector i=l,2 . . . n may have a different form such as 
one having a sine wave form and the other a square wave 
form. However, as already noted, a raised-cosine form is 
well suited given the passenger comfort constraints. 
Consequently, to have independence and still use a raised- 
cosine excitation form, the excitation signal for each redun- 
dant control effector may be applied with different periods. 
This technique can allow two or more redundant control 
effectors to be excited in order to provide for sets of data to 
be gathered from all effectors for a joint optimal solution 
with only one APO algorithm sequence. 

For multiple redundant surfaces, sequential optimization 
could also be used to avoid an identification problem by 
performing the excitation and optimization of one control 
effector first, then set that control effector to its optimal 
position and repeat the procedure for each additional redun- 
dant control effector. Because this will not provide a “joint 
optimal” solution, iteration of this multiple-effectors proce- 
dure should continue for all redundant effectors until con- 
vergence is achieved. The results would then not be achieved 
in what may be considered real time, but near enough to real 
time to be good APO control. 

SECOND 

Data is collected during the excitation of forced response 
maneuver of each one or more redundant control effectors at 
a nominal rate, for example one set of response data per 
second, such that for a raised-cosine excitation of, for 
example, 300 sec., about 300 separate sets of data can be 
acquired, although data acquisition of ten sets per second 
would be preferable for assurance of better optimization. 
THIRD 

For analysis of the sets of response data, the thrust 
estimation is first made for use in equations (1) and (2) 
above in the STATEMENT OF THE INVENTION to cal- 
culate the values of the lift and drag coefficients in those 
equations for each set of data. These sets of data are then 
used to form 300 drag coefficient expansion equations using 
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equation (3), which are then solved for the 5+2(n) unknowns 
in equation (3) (where n=the number of redundant control 
effectors), namely 

Cdo, 4- ? K-i s 6 l0 pt s. 

Solutions can be obtained by regression analysis or other 
comparable techniques. 

10 FOURTH 

Apply the optimal solution (s) to the control efifector(s) in 
a smooth fashion. Use of the first-half of a raised-cosine 
form to reach the optimal value is satisfactory. 

Repeat the optimization algorithm either continuously at 
1 5 fixed flight time intervals or only as needed, i.e., when need 
is determined by changes in flight conditions specified as 
criteria. 

The flow chart of FIG. 4 illustrates the sequence of 
operations for the APO algorithm which recurs under control 
20 of the programmed APO with elapsed time AT or change in 
either altitude AH or Mach number AM. The algorithm 
comprises four steps described above and represented by 
blocks 1 through 4 in the flow chart. Once the optimal h iopt 
is identified and the redundant control effector is positioned 
25 accordingly, the aircraft can then take advantage of the 
optimal (minimum drag) configuration with any number of 
objectives such as but not limited to: 

Minimize Fuel Flow: Trim the aircraft at the same Mach 
and altitude the APO algorithm determined h iopt ; this 
will result in minimum fuel flow. 

Maximize Mach Number: Trim the aircraft at the same 
altitude the APO algorithm determined h iopt and at the 
steady-state fuel flow recorded just prior to the APO 
35 algorithm being requested to determine h iop ^, this will 
result in maximum Mach number. 

Maximize Altitude: Trim the aircraft at the same Mach 
number the APO algorithm determined h iopt and at the 
steady-state fuel flow recorded just prior to the APO 
40 algorithm being requested to determine h iop] ;, this will 
result in maximum altitude. 

Maximum Loiter Time at a Fixed Altitude: Trim the 
aircraft at the same flight condition the APO algorithm 
determined h iopt ; this will result in reduced fuel flow at 
45 that flight condition. Trim the aircraft to a fixed incre- 
ment lower Mach number (e.g., AM -0.02) by reducing 
throttle position (fuel flow). Since this is a new flight 
condition, request the recurrence of the APO algorithm 
to determine a new h iopt which may be imperceptibly 
50 small since the flight condition change was relatively 
small. Repeat these two procedures (reduce Mach num- 
ber and request the recurrence of the APO algorithm for 
a new b iop f ) until fuel flow begins to increase, thus 
locating a h iopt for a minimum Mach number. There 
55 will be that minimum since as Mach is reduced and the 
angle -of- attack increases to maintain lift (in order to 
maintain altitude), a point is reached where drag to 
angle -of- attack is going to be dominant. 

Once the optimal h iopt is identified and the redundant control 
60 effector is moved to the optimal position (minimum drag), 
the APO algorithm (program) can be placed in an idle mode 
until either automatic or crew monitoring of elapsed time or 
change in altitude or Mach number requests a recurrence of 
the APO algorithm. 

65 On an operational flight mission, the APO is turned on by 
the flight crew during any steady-state segment of the 
mission. In any instance of turning on the APO algorithm, 
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the program of the algorithm will first check on whether 
there are any previous optimization results stored. If not, as 
in the case of the APO algorithm being turned on for the first 
time for this particular segment of the mission (e.g., climb, 
cruise, descent or other), the algorithm proceeds to the first 
phase. If the check is positive (YES), the previously stored 
optimal results of the redundant control effector(s) will be 
used as a priori values. This procedure of storing previously 
determined optimal redundant control effector values and 
using the stored values as a starting point when the APO 
algorithm is turned on again will maximize benefits of the 
APO algorithm by utilizing representative optimal values to 
arrive at new optimal values more quickly than starting 
without any a priori values. 

There is a separate 4-dimensional table for each redundant 
control effector(s) value of Mach, altitude, gross weight and 
center of gravity, e.g., as follows: 


DATA TABLE 

Mach 

Altitude 

Gross. Wt. 

c . g . 

min delta 

0.02 

2,000 

10,000 

0.5 

max interp. 

0.10 

10,000 

50,000 

2.5 

Data is stored to 

the nearest min 

delta. The 

maximum 


interpolation allowed to determine an a prioi position is max 
interp. The values in the above table are only representative; 
actualvalues depend on the aircraft in question. 

Note: The redundant control effector(s) can include but is 
not limited to designated effectors of the aircraft such as 
symmetric ailerons, flaps, slats, spoilers, horizontal 
stabilator, elevator, thrust, center of gravity, variable inlet 
and/or engine geometry (for propulsion systems with that 
capability), and thrust vectoring for longitudinal effectors 
and rudder, ailerons, (normal or differential) and differential 
thrust for lateral and directional effectors. 

Although the above process may infer a post-maneuver 
analysis, the analysis can be performed in real time with the 
data acquired to that point in time. In certain applications 
this may provide some benefits. This analysis processing 
then continues for the last part of the raised-cosine excitation 
during which the maneuver is completed and then the entire 
process is repeated until a criteria has been met to assure that 
optimal conditions have been reached (e.g. the criteria that 
the estimated h iopt remains constant over two or more passes 
through the process). 

It would be well to note at this point that the raised-cosine 
excitation form need not begin to immediately diminish in 
amplitude once the maximum amplitude is reached. Instead 
the form may be flat at maximum amplitude for say 200 sec. 
to extend the period of a 300 sec. raised-cosine excitation to 
500 sec. This has the benefit of smoothing the maneuver 
more and allowing more time for the analysis of data 
gathered during the first part of increasing amplitude since 
new sets of data during the flat part need not be processed 
at the times that there is not some significant change from 
the last point of significant change. This is a particularly 
useful benefit as the number n of redundant effectors 
increases beyond a few. 

For normal transport operation in cruise flight, continuous 
optimization is not necessary and in fact would tend to 
increase drag since ideally, once you are at the optimal 
control position of a redundant effector, all other maneuvers 
would be away from that optimal position and therefore 
increase drag. Thus, the optimization algorithm should be 
performed only when either a certain change in flight 
condition criteria is met or a certain time interval has 
elapsed. For both of these, subjective but reasonable criteria 
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would be the following incremental changes: a) AMach= 
0.02; b) Aaltitude=2000 ft.; and Atime=30 min. 

Algorithm solutions can range from continuous to batch 
operation. In either case, all unknown coefficients of equa- 
5 tion (3) are identified including the optimal aileron position, 

a 

®iopt’ 

Values of K 0 and C LminC are assumed to be independent 
of the symmetric aileron position for small deflections. In 
actuality, symmetric aileron deflection would result in a 
10 small change to spanwise lift distribution and, in turn, to 
induced drag characteristics of the wing. This effect would 
cause a small variation in the value of K 0 as a function of 
symmetric aileron position. Sensitivity of the APO algo- 
rithm to this variation has been left for future studies. 

15 The formulation of equation (3) is not unique; the impor- 
tant element being that the first-order effects of aileron- 
induced drag be represented in the C D equation in a plausible 
manner. Care should be taken not to over parameterize the 
problem; independence of the various estimates must be 
20 maintained to provide meaningful results. The actual drag 
reduction is 

ADmg=qSK i (Abf (5) 

Other performance related calculations, such as specific fuel 
25 consumption and range, can be calculated postflight. 
Control Effectors 

A wide range of controls or variables can potentially play 
a role in performance optimization for current and future 
generation aircraft. These controls and variables include 
30 elevator, horizontal stabilizer, outboard aileron, inboard 
aileron, flaps, slats, rudder, e.g., thrust modulation, thrust 
vectoring, and differential thrust (FIG. la) In addition, the 
potential for flightpath control using only differential thrust 
has been demonstrated. Spoilers are probably not an option 
35 for performance optimization because they only increase 
drag. Although if a case exists which requires drag 
modulation, spoilers are a viable controller. 

Note that delta-wing aircraft configurations generally 
have less optimization potential than tail-configured aircraft 
40 because delta wings have fewer independent control effec- 
tors. Fewer effectors reduce the potential for optimization. 
The main difference is that there is no independent horizon- 
tal stabilizer-elevator for delta-wing configurations, thus 
removing a major potential for wing optimization. These 
45 differences do not imply that more sophisticated wing 
leading- and trailing-edge (TE) devices, which would permit 
some degree of camber optimization, could not be imple- 
mented. A canard can significantly increase the optimization 
potential for conventional and delta-wing configurations. 

50 Instrumentation 

Successfully implementing a performance optimization 
algorithm requires high-quality, sensitive instrumentation. 
Fortunately, the instrumentation being implemented in 
today's most advanced transports for FMS operation 
55 (trajectory and navigation control) is very good. Although a 
large number of cost functions or variables exists that could 
conceivably be used for optimization, only a few basic 
aircraft measurements are required for cruise drag minimi- 
zation. To minimize fuel flow at constant Mach number and 
60 altitude conditions requires accurate fuel-flow indications, 
such as either fuel flow, fuel valve position, throttle position, 
or thrust. 

In-lieu of direct fuel-flow measurements, engine pressure 
ratio (EPR) measurements combined with a representative 
65 engine model, which is a function of flight condition, will 
provide accurate incremental fuel-flow or thrust results. 
Although in absolute terms the accuracy required would be 
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demanding, the optimization problem only places demands 
on perturbation accuracy, which is not affected by biases. 
Maximizing velocity for constant altitude and fuel flow 
requires accurate perturbation measurements of velocity, 
flightpath acceleration, or both. 

Real-Time Drag Minimization 

Preferably, performance optimization could be accom- 
plished using responses to pilot or autopilot and FMS 
commands. However, with tight pitch-rate, pitch-attitude, 
and altitude- and Mach-hold control laws, external 
environment-based disturbances and associated responses 
would, generally, be small. As a result, forced response of 
redundant controls by unique excitation is required to ensure 
identifiability. The requirement for forced excitation must be 
tempered by the additional requirement that neither handling 
nor ride qualities are noticeably impacted. In turn, this 
requirement dictates the range of excitation periods and 
amplitudes. 

Aircraft Simulation Model 

The APO algorithm was evaluated using a simulation of 
a first-generation wide -body transport. The simulation is 
high fidelity and covers the full aircraft envelope. The 
primary control system has simple rate feedback, while the 
altitude- and Mach-hold autopilot modes were designed with 
feedback to the stabilator and throttle, respectively. The 
aileron drag characteristics were modified (based on simu- 
lation and flight data of similar configurations) to provide a 
quadratic drag variation with symmetric aileron displace- 
ment. All the simulation runs were initiated at Mach 0.83, at 
an altitude of 37,000 ft, and at a weight of 408,000 lb. Fuel 
burn was simulated as a function of thrust. Low-order thrust 
versus throttle-lever dynamics are in the simulation. Light 
turbulence was used for all runs to provide a realistic 
signal- to -noise ratio. FIG. 3 shows the variation of thrust 
required for trimmed steady-state flight as a function of 
symmetric aileron displacement. The maximum drag reduc- 
tion of 364 lb (1.5 percent) occurs at approximately 4.5° 
symmetric aileron deflection. 

Results and Discussion 

FIG. 5 presents the forced excitation response, with 
altitude- and Mach-hold modes on, to a raised cosine (1.0- 
cos (cot)) symmetric aileron command with a 300-sec dura- 
tion. (The co is the angular frequency.) The simulated 
responses used in the analysis are flightpath axes 
accelerations, angle of attack, thrust, weight, and symmetric 
aileron deflection. For this analysis, it was assumed that the 
angle of attack is calculated from pitch attitude INS veloci- 
ties and true airspeed. In addition, the flightpath accelera- 
tions are calculated from body-axes accelerations rotated 
through angle of attack, and thrust is obtained as a function 
of EPR and flight condition. The variation of C D with 6, does 
not present a clear picture of the minimum drag point (FIG. 
6 ). 

FIG. 7 presents &C D =C D -K 0 {C L -C LminCo } 2 with 5,. being 
the sole symmetric aileron, the sole redundant control effec- 
tor. Correcting for the C L variations produces a much clearer 
picture of the drag minimization process. The quadratic 
variation of AC^, is clear, but there still is a significant 
difference between increasing and decreasing 6,- commands. 
(Note that the C D and AC^, incremental scales are the same 
for FIGS. 6-9 for ease of comparison.) 

A close look at the time histories reveals that Mach 
number variations exist which could, in turn, contribute a 
change in drag as a function of Mach number via a C Dm 
effect. To improve the accuracy of the analysis, a linear 
C^AM term was added to Equation 3 with Chadded to the 
list of variables to be estimated. FIG. 8 presents the AC^ 


variation (including the Mach number effect correction) with 
5-=5 fl (for symmetric aileron) and the quadratic variation is 
now well-defined and agrees with the minimum drag point 
of FIG. 4. 

5 The variation of drag with the b iopt variation can be 
modeled as a quadratic function for the purposes of locating 
the local minimum of even a potentially complex function. 
It is not required that the true drag associated with the h iopt 
variation be quadratic but only that a quadratic representa- 
10 tion of the local minima be reasonably representative of the 
true minima. This modeling will assure that the local minima 
is found. 

Sensitivity to Baseline Data 

The good results presented in FIG. 8 may not be all that 
15 surprising in view of the fact that the simulation output 
variables are perfect with, at most, some white noise effects 
manifested from the turbulence. Assurance of analysis 
insensitivity to all known effects must be verified. 

To minimize the variation or inaccuracy of the estimates, 
20 or both, a priori values of K 0 and C Lmin c D are used. These 
parameters characterize the quadratic nature of the C L /C D 
variation. Separate parametric variations of these a priori 
parameters of ±25 percent and ±50 percent on K 0 and 
C L minc D respectively were performed. These variations pro- 
25 duced less than 0.1° variation on the optimal symmetric 
aileron value (4.5°). As a result, the optimal solution appears 
insensitive to these a priori parameters. 

Measurement Bias Effects 

For absolute performance analysis, measurement bias is a 
30 limiting factor on analysis accuracy. However, the formu- 
lation of this APO analysis algorithm is designed to be 
insensitive to measurement bias. The following biases 
(applied one at a time) produced less than 0.1° variation on 
the optimal symmetric aileron value: 

35 


Parameter 

Bias 

Angle of Attack 

0.5° 

Acceleration along the flightpath 

0.02 g 

Acceleration (upward) normal to the flightpath 

0.02 g 

Net aircraft thrust 

1500 lb 

Aircraft gross weight 

10,000 lb 


A bias on the symmetric aileron measurement produces an 
45 equivalent change on the optimal solution, but in true terms, 
the solution is unaffected. Bias insensitivity is significant 
relative to the analysis because measurement biases are 
common. Note also that these findings apply in spite of the 
fact that the equations are nonlinear. 

50 Resolution Effects 

Operating near or at the resolution limits of instrumenta- 
tion is potentially a serious problem, and the analysis 
procedure must be insensitive to these quantization effects. 
Formulation of this APO analysis algorithm is not overly 
55 sensitive to resolution resulting, in part, from the regression 
technique employed. The results of FIG. 8 are repeated in 
FIG. 9 but with the following resolution set: 


Parameter 

Bias 

Aileron position 

0.10 

Angle of Attack 

0.10 

Acceleration along the flightpath 

0.002 g 

Acceleration (upward) normal to the flightpath 

0.002 g 

Net aircraft thrust 

150 lb 

Mach number 

0.001 lb 
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The analysis with resolution produced less than 0.2° varia- 
tion from the optimal symmetric aileron value case of FIG. 

8 . 

The algorithm appears to be most sensitive to thrust 
resolution. An increase in the thrust resolution to 300 lb 5 
changes the optimal symmetric aileron value by 0.4°. In 
general, however, the optimal solution is relatively insensi- 
tive to these resolution effects. 

Thrust Accuracy 

Although thrust is not measured, it is determined based on 
EPR and flight condition measurements and a representative 
engine model. The thrust calculations will tend to be the 
least accurate of all inputs to the analysis process. Constant 
errors in the thrust level are not a problem. However, thrust 
is based on interpolation of steady-state thrust tables; hence, 
inaccuracies caused by lack of modeled engine dynamics, 15 
whether they accrue from throttle lever motion or atmo- 
spheric effects, will occur. A cursory evaluation of this effect 
was conducted by calculating a thrust value to be fed into the 
analysis as a linear variation with throttle lever. FIG. 10 
shows the variation of thrust with throttle lever for the data 20 
of FIG. 5 along with a linear fit of the data. Using this 
relationship in the analysis produced an error of less than 
0.1° in the optimal symmetric aileron solution. 

Minimizing fuel flow at constant Mach number and 
altitude conditions requires accurate fuel-flow indications, 25 
such as fuel flow, fuel valve position, throttle position, or 
thrust. In-lieu of direct fuel-flow measurements, EPR mea- 
surements combined with a representative engine model, 
which is a function of flight condition, will provide suffi- 
ciently accurate incremental fuel-flow or thrust estimates. 30 

A more realistic and interesting situation is the sensitivity 
of the solution to an error in the slope of the thrust as a 
function of the throttle of the linear fit (FIG. 10 ). A bias of 
the linear relationship is addressed above under the Mea- 
surement Bias Effects heading. An error in the slope (in 35 
addition to a bias error) would be representative of a 
miss-modeling of the engine characteristics via a table 
look-up process. FIG. 10 also shows a linear thrust-throttle 
lever variation with a 10-percent increase in slope and with 
the constant adjusted so that the linear relationship has the 40 
correct thrust level at the trim flight condition. The constant 
has the same effect as a bias error and, therefore, is a 
reasonable assumption. The 10-percent error in slope pro- 
duced a negligible (less than 0.1°) difference in the optimal 
aileron solution. A 20-percent error produces a 0.4° differ- 45 
ence in the optimal aileron solution. In referring back to FIG. 

5, even a 0.5° error in the optimal solution is insignificant 
because the variation of drag with aileron is shallow. 
Hardware Implementation 

Selected aircraft have the hardware (symmetric aileron 50 
trim deflection capability) required to perform onboard 
performance optimization and, therefore, would only require 
a relatively simple set of optimization software. Algorithm 
redundancy is not required because the algorithm is a 
nonsafety-of-flight system. This nonsafety-of-flight aspect 55 
can be assured by having APO in a discretionary mode with 
very limited rate and position authority. The algorithm can 
be a completely independent set of code and, therefore, 
avoid the issues of integration with the FMS. The algorithm 
can be thought of as a slow, limited authority trimmer. 60 

For aircraft that do not currently have the hardware 
capability of moving the outboard ailerons symmetrically, a 
relatively simple modification consisting of adding a low 
frequency, limited authority, trim actuator in series with the 
mechanical command to the outboard aileron actuator would 65 
suffice. The slow actuator rate plus limited authority mini- 
mizes safety-related issues. 


Related Issues 

The forced excitation requirement of real-time adaptive 
optimization generally attracts concern and, therefore, some 
discussion is in order. For the very steady conditions of 
cruise flight optimization, forced response excitation is the 
only means of performing identification and adaptive opti- 
mization. No other means of identifying the desired char- 
acteristics of the aircraft exist. The low frequency raised 
cosine excitation was intended to minimize interaction with 
the autopilot inner-loops and to minimize the effect on ride 
qualities. This function results in a negligible incremental 
normal acceleration (*0.001 g) and precludes concern for 
other aircraft wear and tear, such as control surfaces and 
engines. Although minimal interactions between the explicit 
excitation and the inner loop autopilot modes were 
observed, other aircraft, control surface combinations, or 
both, could produce degraded inner-loop performance. In 
such instances, inner-loop controller lead as a function of the 
explicit excitation function could be used to minimize 
inner-loop performance degradation. 

Because aircraft specific variations play a significant role 
in the actual amount of performance improvement 
accruable, using previous optimality results as initial con- 
ditions can speed up optimality convergence for subsequent 
flights. 

Different flight conditions have also been evaluated, and 
it appears that very little, if any, algorithm tuning is required. 
Exploratory algorithms have also been evaluated for opti- 
mization in the climb- to -altitude flight segment. 
Concluding Remarks 

This invention described presents a real-time, adaptive- 
optimal performance algorithm for application to subsonic 
transports. Preliminary simulation results indicate the 
approach is very promising. The algorithm implementation 
is simple and appears to have robust performance charac- 
teristics. 

Because the performance optimization problem is search- 
ing for small benefits, instrumentation may appear to be a 
critical factor. However, realistic instrumentation and mea- 
surement effects have been evaluated, demonstrating that the 
algorithm is insensitive to these effects. 

While a simulation model demonstrates that the algorithm 
looks promising, real world effects cannot be predicted or 
simulated accurately. Nevertheless, it can be predicted and 
simulated that application of APO to any fly -by-wire aircraft 
which already has electronic control over desired redundant 
control effectors, requires only software additions which 
would logically be placed in the FMS, and that the pilot 
would interface the APO software through the Control and 
Display Unit (CDU) located in the cockpit. To that extent a 
simulation model demonstrates that APO is applicable to 
any fly-by-wire aircraft. Aircraft which do not have elec- 
tronic control over the redundant control effectors, need only 
to be modified by installation of electrically commanded 
actuators (or equivalent) somewhere within the system con- 
trolling the particular redundant control effectors of interest. 
Of course, this would also require the software modifica- 
tions in the FMS discussed above for fly-by-wire aircraft. A 
flight evaluation of the adaptive performance optimization 
algorithm is underway to research the issues of operational 
use; benefits can only be determined in-flight. However, it 
can now be stated that the application to supersonic trans- 
ports will have a significantly greater benefit than for 
subsonic transports, because the supersonic transport must 
cover a much larger envelope, has a variable geometry 
propulsion system, and thus provides more opportunities for 
optimization. 
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Although particular embodiments of the invention have 
been described and illustrated herein, it is recognized that 
modifications may readily occur to those skilled in the art. 
Consequently, it is intended that the claims be interpreted to 
cover such modifications and equivalents thereof. 5 

What is claimed is: 

1. A method for providing in-flight adaptive performance 
optimization of an aircraft in a steady climb, cruise, turn or 
descent flight condition that uses autopilot and navigation 
modes in a normal manner or by the pilot’s direct control, 10 
where said optimization is based on real-time measurements 
and calculations of incremental drag from a forced response 
maneuver of one or more redundant control effectors that 
affect aircraft drag defined as those control effectors that are 

in excess of the minimum set of control effectors that affect 15 
aircraft drag required to maintain said steady flight condition 
in progress, comprising the steps of first applying excitation 
to said one or more redundant control effectors simulta- 
neously in a smooth raised-cosine form over a finite interval, 
concurrently gathering sets of measurements made during 20 
said excitation, next calculating lift and drag coefficients C L 
and C D from two equations, one for each coefficient, and 
then calculating an optimum position b iopt of said one or 
more redundant control effectors from a third equation 
which is an expansion of C D as a function of parasitic drag, 25 
induced drag, Mach, altitude drag effects, and control effec- 
tor drag, said third equation comprising a quadratic variation 
of drag with positions b i of redundant control effectors i=l 
to n, solving said third equation for b iopt in a separate step, 
and finally using the value b iopt to set the position of each 30 
control effector i for optimum performance of said aircraft 
throughout the remainder of said steady flight condition, 
unless monitored flight conditions change by some prede- 
termined minimum amount, or a predetermined minimum 
flight time lapses, or flight conditions have changed an 35 
amount perceptible by a flight crew monitoring instruments 
in the cockpit of said aircraft, in response to any one of 
which occurs, said method for providing in-flight adaptive 
performance optimization is repeated. 

2. A method as defined in claim 1 for providing in-flight 40 
adaptive performance optimization of an aircraft after said 
steady flight condition has been reached wherein, 

while applying forced response excitation to said one or 
more redundant control effectors simultaneously, col- 
lecting sets of data during the forced response 45 
excitation, performing optimization analysis to deter- 
mine lift and drag coefficients C L and C D for each set 
of data by solving the following equations: 

C L =Lift/( ? 5)={ m z//7 -2r, sin (a-ri)}/(?S) 50 

C 0 =Drag/(#S)={27;. cos (a-r\)-WA xfp }!(qS) 

where: q=dynamic pressure, lb/ft 2 , 

S=wing reference area, ft 2 , 

A=angle of attack, radians, 

W=aircraft gross weight, lbs, at the start of said method 
A zf =acceleration normal to flight path, (positive 
upward), 

A^y =aircraft acceleration along the flight path (positive 60 
iferward), 

r|= inclination of engine thrust relative to aircraft fuselage, 
radians, 

Z=summation of the bracketed quantity over the number 
of engines of the aircraft from i=l to N, solving for 65 
estimates of b iopt for each redundant control effector i 
by calculating unknown variables 


Cj) 0 , K 0 , CijninCj^Cj)^ 4- C Dfjp , K; S 

in the following equation: 

Q> = C Do 4- K 0 \C L - CiminCj) } 2 + C Dm AM 4- C DHp A H p 4- XKiVi - 6 iopt f 

where C^ 0 =minimum drag coefficient, 

C Lminc D = C L at minimum C D , 

(^^coefficient of drag caused by Mach number, 

C ' D H p = coefficient of drag caused by altitude, 

AM=change in Mach number from the time of turning the 
APO on, 

AHp=change in altitude from the time of turning the APO 
on, 

K 0 and K, are drag equation coefficients, 

5 —position of redundant control effector i, 

5 IO/3r =opti ma l position of redundant control effector i, 
2=summation over a number n of redundant control 
effectors from i=l to n and finally solving for b iopr 

3. A method as defined in claim 2 wherein the equation for 
is a quadratic variation of C^, with 8,-. 

4. A method as defined in claim 2 wherein said one or 
more redundant control effectors of said aircraft comprise 
symmetric ailerons, flaps, slats, spoilers, stabilizer, elevator, 
thrust, center-of -gravity and thrust vectoring, all for longi- 
tudinal effectors, and rudder, normal or differential ailerons 
and differential thrust, all for lateral or directional effectors. 

5. A method as defined in claim 2, wherein to reduce the 
number of unknowns, a priori values of K 0 and C LrninCo are 
utilized to solve for optimal b iopt estimates defined by the 
equation for C^, from which, at the end of the analysis of all 
sets of data the most optimal estimate b iopt is determined for 
the position of redundant control effector i until said method 
is repeated for an update of the setting of redundant effector 
positions to maintain optimal adaptive performance opera- 
tion of the aircraft. 

6. A method for providing in-flight adaptive performance 
optimization of an aircraft as defined in claim 1 wherein the 
number of redundant control effectors is equal to 1, said 
method comprising using a forced response maneuver of 
said redundant control effector to collect aircraft response 
data to forced maneuvers of limited extent to search for the 
optimal redundant control effector position, b iopt that mini- 
mizes drag by solving the following lift and drag equations: 

C L =lMl{qS)={WA zf -YT i sin(a-ri)}/(?S) 

C^Drag/foSHZT, cos(a-^-(W)A xfp }/(qS) 

where: q=dynamic pressure, lb/ft 2 , 

S=wing reference area, ft 2 , 

A= angle of attack, radians, 

W=aircraft gross weight, lbs, at the start of said method 
A zf = acceleration normal to flight path, (positive 
upward), 

.A^j^y = aircraft acceleration along the flight path (positive 
&rward), 
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T| = inclination of engine thrust relative to aircraft fuselage, 
radians, 

Z=summation of the bracketed quantity over the number 
of engines of the aircraft from i=l to N, solving for an 
estimate of h iopt for a redundant control effector by 
calculating the other unknown variables 

Cdoi K-oi CfjninCpCDm + Cj)^^ , 

in the following equation: 

Cd = C Do + K 0 \C l - C Lmin c D } 2 + C Dm AM + C Dh ? LH p + K\ (t5 x - 6 lopt f 

where ^Do =minimum drag coefficient, 

C Lminc= C L at minimum C D , 
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C Z3w =coefficient of drag caused by Mach number, 

Co Hp = coefficient of drag caused by altitude, 

5 

AM=change in Mach number from the time of turning the 
APO on, 

AH p =change in altitude from the time of turning the APO 
on, 

io K 0 and K a are drag equation coefficients, 

S -^position of said redundant control effector, 

6 lc ^ r =opti ma l position of said redundant control effector, 
said aircraft response data collected for analysis using the 
above equations being acquired at a steady flight condition 
15 comprising relatively constant Mach number, altitude and 
heading accomplished by any combination of autopilot and 
pilot control. 



